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Abstract 

The refraction of linearly polarized plane waves into a half-space occupied by a material moving 
at constant velocity was studied by directly implementing the Lorentz transformations of electric and 
magnetic fields. From the perspective of a co-moving observer, the moving material was a spatially 
local, pseudochiral omega material. Numerical studies revealed that whether or not negative refraction 
occurrs in the moving material depends upon the speed of movement as well as the angle of incidence 
and the polarization state of the incident plane wave. Furthermore, the phenomenons of negative phase 
velocity and counterposition in the moving material were similarly found not to be Lorentz covariant; 
both phenomenons were also found to be sensitive to the angle of incidence and the polarization state of 
the incident plane wave. 

PACS numbers: 03.30.+p, 03.50.De, 41.20.Jb 

1 Introduction 

The behaviour of plane waves at a planar interface between two disparate homogeneous mediums is a central 
topic in both fundamental and applied electrodynamics. In particular, the phenomenon of negative refraction 
[1, 2, 3] has been the subject of intense research efforts for the past ten years, following experimental reports of 
this phenomenon in certain metamaterials [4, 5]. Much of this effort has been motivated by the development 
of novel metamaterials, but negative refraction also arises in certain minerals [6] and in biological structures 
[7] . In addition, the prospects of negative refraction arising in rclativistic scenarios — such as in uniformly- 
moving materials [8, 9] or in strong gravitational fields [10, 11, 12, 13] — is a matter of astrophysical and 
astronomical significance. 

In the following we consider negative refraction induced by uniform motion. Earlier work relating to 
this topic relied on the Minkowski constitutive relations to describe the moving medium in a nonco-moving 
inertial reference frame [14, 15, 16], per the standard textbook approach [17]. However, this approach is only 
appropriate to materials which are both spatially and temporally local [18]. More recently, studies based 
on the uniform motion of realistic materials have been undertaken, using an approach in which the Lorentz 
transformations of the electric and magnetic fields arc directly implemented [8, 9]. These studies revealed 
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that the phenomenons of negative phase velocity and counterposition 3 — which are closely allied to negative 
refraction and similarly associated with certain metamaterials and relativistic scenarios — are not Lorentz 
covariant. Here we address the hitherto outstanding question: is negative (or positive) refraction Lorentz 
covariant? By means of a numerical analysis based on a uniformly moving pscudochiral omega material, we 
demonstrate that the answer to this question is 'no'. 

In the notation we adopt, 3-vectors are in boldface with the " symbol denoting a unit vector. Double 
underlining signifies a 3x3 dyadic (i.e., a second rank Cartesian tensor) and the identity 3x3 dyadic is 
written as/ = xx + yy + zz. The operators Re and Im deliver the real and imaginary parts of complex 
quantities; and i = The permittivity and permeability of free space are e and fx , respectively, with 

Co = l/^/CoMo being the speed of light in free space. 



2 Refraction into a moving pseudochiral omega material 
2.1 Pane wave analysis 

Our attention is focussed on a spatially local, homogeneous material, characterized by the frequency-domain 
constitutive relations 

D' = i • E' + f ' • H' 1 



B' = C' • E' + fj • H' 

in the inertial reference frame Herein, the 3x3 constitutive dyadics 
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This is a bianisotropic, Lorentz-reciprocal [19] material, known as a pseudochiral omega material [20]. Con- 
stitutive relations of this form have been used to describe certain negatively refracting metamaterials, as- 
sembled from layers of split-ring resonators [21]. Several different designs of metamaterials are based on this 
general configuration [22, 23, 24]. As the pscudochiral omega material is presumed to be dissipative, the 
constitutive parameters e' x y 2 , f and fj,' x y z are complex- valued functions of the angular frequency u'. 

Suppose that the pseudochiral omega material fills the half-space z > 0, while the half-space z < is 
vacuous. The inertial reference frame £' translates at constant velocity v = vv with respect to the inertial 
reference frame E, in the plane of the interface z — 0. In keeping with an earlier study [9], we take v = x. 
The Lorentz transformations [17] 



E 
B 
H 
D 



(E' • v) v + 7 [{I - vv) • E' — v X B'] 

v x E' 



with the real valued scalars 



(B' • v)v + 7 
(H' • v) v + 7 [(I- vv) 
(D' • v) v + 7 (1 - vv) 

1 



B' 
H' 
•D' 



v x D'] 
v x H' 



(3) 



/3 



v 

Cn 



(4) 



'See Sec. 2.2 for definitions of negative refraction, negative phase velocity and counterposition. 



2 



relate the electromagnetic field phasors in the frame E to those in the frame E'. 

Now suppose that the vacuous half-space z < contains a line source, which is stationary with respect 
to the frame E. The source extends infinitely in directions parallel to the y axis, and it is located at a great 
distance from the interface z = 0. Let us consider one plane wave incident on the interface z = 0, as a 
representative of the angular spectrum of plane waves launched by the source. With respect to the frame E, 
this plane wave is described by the electric and magnetic field phasors 



Ej = e 4 exp [i (k; • r - cut)] 
Hi = hi exp [i (kj • r — cot)] 



z<0. 



Herein, the wavevector 
with the real-valued scalar 



(tx 



k cos( 



k = k a sin# e (— k , k ) , 



the free-space wavenumber k„ = ujy/e fi and to being the angular frequency with respect to E. 
With respect to the frame E', the incident plane wave is represented by 

E<=e£exp[i(k^r'-a/i')] 1 
H^exp^.r'-c^')] J ' 

wherein the phasor amplitudes {ej, h'j} are related to {e^h^} via the Lorentz transformations (3), while 



(5) 

(6) 
(<) 

(8) 
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k • v + — 



I- vv 



k' 



r = [L+ (7- 1) *v] 'r' + jvt' 
lo = 7 (u)' + k'i • v) 
( 1 v • r' 

* = 7 K + — 5- 



(9) 



The incident plane wave gives rise to two refracted plane waves in the half-space z > 0, and one reflected 
plane wave in the half-space z < 0. In the frame E', the refracted plane waves are represented by the electric 
and magnetic phasors 



E' t = e' tj exp [i (k' tJ ■ r> - w'f )] 



h: 



h',- exp [i (k' f • r' 



to't' 



z > 0, 



(i = i,2), 



wherein the wavevectors 



k 'tj = • x) x+ z, (j = 1,2) 



(10) 



(11) 



comply with Snel's law [17]. The wavevector components k' Z j, as well as the relationships between the phasor 
amplitudes ejy and h^-, are deduced by combining the constitutive relations (1) with the source-free Maxwell 
curl postulates in E' [17]. We find [29] 



k' zX = uj'^e fi a 
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Notice that since k' z j are generally complex-valued, the refracted plane waves are nonuniform. 
The reflected plane wave is represented by the electric and magnetic phasors 



E; = e' r exp [i (k; • r' - w't')] 

h; = h; exp [i (k; ■ r' - w't')] 



z<0, (13) 



in the frame £', with the wavevector of the reflected plane wave being 

k' r = (kj >x)x-k' zr z. (14) 

The source-free Maxwell curl postulates in £' yield an expression for k' zr and relationships between the 
phasor amplitudes e' r and h' r . 

By invoking the standard boundary conditions across the plane z — 0, i.e., [17] 



fa + e' r ) ■ x = e' tJ ■ x (e' t + e' r ) ■ y = e' tj ■ y 
(h' ( + h;) • x = h' tj ■ x (h^+h;).y = hj..y 



z = 0, (i = 1,2), (15) 



the phasor amplitudes {e£.,hj,} and {ej-,hj-} can be found. The reflected and the refracted plane waves in 
the frame E may then be deduced by applying the Lorentz transformations (3) and (9). 

We represent the refracted plane wave in the frame £ by the electric and magnetic field phasors 



E t = e t j exp [i {k tj ■ r - tot)] 
H t = h tj exp [i (k tj ■ r - ut)] 



z>0, (7 = 1,2), (16) 



wherein the wavevectors 

k tj =KX + k zj z, (J = 1,2), (17) 
while the corresponding reflected plane wave is represented by 



E r = e r exp [i (k r • r — u>t)] 
H r = h r exp [i (k r • r — ut)] 



(18) 



with the wavevector 

k r = Kx- k cos 9 z. (19) 
2.2 Negative refraction, negative phase velocity and counterposition 

We highlight three phenomenons which the uniformly moving pseudochiral omega medium can support: 
namely, negative refraction, negative phase velocity and counterposition. All three phenomenons are closely 
associated with certain mctamaterials [29]. Wc consider these phenomenons referred to the frame S. 

First, in the scenario under consideration here, negative refraction occurs when the real part of the 
wavevector ky casts a negative projection onto the positive axis z; i.e., when Rek Z j < [2]. Second, the 
phase velocity of the refracted wave is classified as negative when the real part of the wavevector k t j casts 
a negative projection onto the corresponding cycle-averaged Poynting vector Pj [1]. For our purposes here, 
it suffices to consider Pj in the limit |r| — > 0, as provided by 

Pjl, rM ) = ^[(Ree tj )x (Rehy) + (Ime y )x (Imhy)], (i = 1,2). (20) 

The signs of the square roots in eqs. (12) are selected in order to ensure that Pj casts a positive projection onto 
the positive z axis. Thereby, energy flow directed into the z > half-space is assured. Third, counterposition 
occurs when the real part of the wavevector k t j and the corresponding cycle-averaged Poynting vector Pj 
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are oriented on opposite sides to the normal to the interface z = [25, 26, 27]. Because here Pj • z > 0, the 
conditions for counterposition are 

(Rek tJ ;• z) (Pj • x) < for/t>0 
(Rek tJ • z) (Pj • x) > for « < 

For uniform planewavc propagation in an isotropic dielectric material, negative refraction and negative 
phase velocity are effectively synonymous [1, 2]. However, this is not the case for more complex materials 
[28]. In particular, for the pseudochiral omega material considered here, it has been established that negative 
refraction, negative phase velocity and counterposition are independent phcnomcnons [29]. 

2.3 Numerical studies 

We now explore the phenomenons of negative refraction, negative phase velocity and counterposition for a 
specific numerical example, chosen to allow direct comparison with an earlier study [29] . Let the pseudochiral 
omega medium occupying the half-space z > be characterized by the constitutive parameters: e' x = 
0.1 + 0.031 e' y = 0.14 + 0.02i, e' z = 0.13 + 0.07i; fx' x = -0.29 + 0.09z, fj,' y = -0.18 + 0.031 fj,' x = -0.17 + 0.6i; 
and £ = 0.11 + 0.051 

In the following we consider two polarization states for the incident plane wave. As described in the 
Appendix, the incident s-polarization state, as characterized by 

ei = a s s = a s y, (22) 

gives rise to the refracted plane wave with wavevector k t i, whereas the incident p-polarization state, as 
characterized by 

e.; = a p p + = a p (— cos 9 x + sin 9 z) , (23) 

gives rise to the refracted plane wave with wavevector kt2 ■ The corresponding reflectances and transmittances 
are presented in the Appendix. 

The orientation angle of the real part of the wavevector k tJ -, (j = 1, 2), as defined by tan -1 (k/Rc k z j), is 
plotted in Fig. 1 versus the relative speed (3 £ (—1, 1) for the angle of incidence 9 G {0°, 5°, 25°}. For normal 
incidence, we see that Re k tJ - is also normal to the z = interface, for both incident s- and p-polarization 
states. However, for 9 = 5° and 9 = 25°, the orientation of Rek t j relative to the interface normal is highly 
sensitive to j3, for both incident s- and p-polarization states. For example, for 9 = 5° with e.j = s, negative 
refraction occurs for — 1 < (3 < —0.03 and for 0.21 < j3 < 1, and the refraction is positive otherwise. 
Similarly, for 9 = 25° with e» = s, and for 9 £ {5°, 25°} with = p + , whether the refraction is negative or 
positive depends upon the relative speed (3. The real part of k tJ - is oriented normally to the z = interface 
in the limits (3 — > ±1, for both j = 1 and 2. 

In Fig. 2, the quantity Pj | . ■ ^ • Reky, (j = 1,2), which determines whether the phase velocity is 

positive or negative, is plotted against relative speed (3 € (—1, 1) for the angle of incidence 9 £ {0°, 5°, 25°}. 
The numerical results echo those of Fig. 1 insofar as the sign of the phase velocity is highly sensitive to (3 
for both states of incident polarization, for 9 £ {5°, 25°}. Unlike the case of negative/positive refraction, 
the sign of the phase velocity is highly sensitive to (3 even at normal incidence, for both the s- and the 

p-polarization states of the incident plane wave. In the limits (3 — > ±1, Pj | ■ r ■ ^ becomes oriented parallel 

to the interface z = for both the s- and the p-polarization states of the incident plane wave and the phase 
velocity of the refracted plane wave therefore becomes orthogonal. 

Lastly, we turn to the quantity (Pj • x) Re (k f j • z), which determines whether or not the cycle-averaged 
Poynting vector and the real part of the wavevector for the refracted plane wave are counterposed, per the 
conditions (21). This quantity is plotted against relative speed (3 £ ( — 1,1) for the angle of incidence 
9 £ {0°, 5°, 25°} in Fig. 3. In a manner similar to that represented by Figs. 1 and 2, we see that the sign of 
(Pj • x) Re (k t j • z) is highly sensitive to f3 for both j = 1 and 2. Furthermore, for all angles of incidence 
considered with e,; = s and with e, = p + , counterposition arises in the limit /3 — > 1 whereas it does not arise 
in the limit j3 — ► — 1. 



(i = i,2). 



(21) 



5 



3 Concluding remarks 



The main conclusions to be drawn from this study are that, for the uniformly moving pseudochiral omega 
material under consideration, 

(i) negative refraction is not Lorentz covariant; 

(ii) negative phase velocity is not Lorentz covariant; and 

(iii) counterposition is not Lorentz covariant. 

While conclusions (ii) and (iii) are consistent with findings previously reported for a uniformly moving 
isotropic dielectric material [9, 14], conclusion (i) reveals a previously unknown facet of refraction. This has 
far reaching consequences for researchers exploring the refraction of electromagnetic waves in astrophysical 
applications, as well as those investigating fundamental aspects of electromagnetic theory. In particular, the 
absence of Lorentz covariance of negative refraction further vindicates the 3+1 approach to electromagnetic 
wave propagation in strong gravitational fields, which has recently been usefully exploited to elucidate 
gravitationally induced negative phase velocity in vacuum [10, 11, 12, 13]. 



Appendix 

Here we present the reflectances and transmittances for the reflection-transmission scenario described in 
Sec. 2. In terms of linearly polarized states, let the amplitude of the electric phasor of the incident plane 
wave be expressed as 

e, = a s s + a p p + , (24) 

where the unit vectors s = y and p + = — xcos(9 + zsin#. The corresponding vector amplitude for the 
reflected plane wave may be written as 

e r = r s s + r p p-, (25) 
where the unit vector p„ = xcos# + zsin^, while for the refracted plane wave we have 

e i = t 1 s + t 2 p, (26) 

where the vector p lies in the plane of incidence and satisfies p • p* = 1. The relative amplitudes r s-p /a s-p 
and ti^/a s _ p may be deduced by following the strategy described in Sec. 2.1; i.e., by transforming to the 
frame then invoking the boundary conditions (15), and finally transforming back to the frame S. By so 
doing, we find that 

a s = =^> r s = ti = I 

(27) 

a p = => r p = t 2 = J 

Thus, there are no cross-polarization terms. 

For the particular example considered in Sec. 2.3, the reflectances \r s /a p \ 2 and \r p /a p \ 2 , and the transmit- 
tances |ti/a s | 2 and \t2/a p \ 2 , are plotted against relative speed (3 G ( — 1, 1) in Fig. 4 for the angle of incidence 
9 e {0°, 5°, 25°}. We note that, for all angles of incidence considered, the reflectance \r p /a p \ 2 tends to unity 
in the limits (3 — > ±1 whereas this is not the case for the reflectance \r s /a s \ 2 . In a similar vein, for all 
angles of incidence considered, the transmittance \t\/a s \ 2 tends to zero in the limits (3 — ► ±1 whereas the 
transmittance |^2 / a p | 2 does not. 
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Figure 1: The orientation angle of Rek tJ (in degree) plotted against relative speed (3 <E (—1,1) for the angles 
of incidence 8 = 0° (broken curve, red), 6* = 5° (solid curve, green), and 8 = 25° (broken dashed curve, blue). 
Plots are shown for the s-polarization state (j = 1) and the p-polarization state (j = 2) of the incident 
plane wave. 




Figure 2: As Fig. 1 except that the quantity plotted against is Pj • Rek tJ - (normalized), (j = 1,2). 
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Figure 3: As Fig. 1 except that the quantity plotted against (3 is (Pj • x) Re (ky • z) (normalized), (j = 1, 2). 
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